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® Verfahren zur in-vivo Bestimmung der Compliance-Funktion und des systemischen Blutflusses eines 
Lebewesens und Vprrichtung zur Durchfuhrung der verfahren 

@ Zur individuation in-vivo Bestimmung der Compliance- 
Funktion C(p) « dV/dp des GefaSsystems nach einer Herz- 
kammer eines Lebewesens aus dem Blutdruck pit) und ei- 
nem Referenz-Herzzeitvolumen HZVref, wird 

a) der Druck p(t) in der Aorta oder in der Nahe der Aorta 
laufend bestimmt; 

b) ein Referenzwert fur das Herzzeitvolumen HZVref ver- 
wendat; 

c) aus dem Blutdruck p(t) der mittlere Blutdruck MAD be- 
rechnet; 

d) der systemischa Widerstand R des Korpers berechnet 
als 
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wobei ZVD ein beliebiger zentralvenoser Druck ist, der er- 
mittelt Oder geschatzt wird; 

e) aus dem Blutdruck mindestens die erste zeitiiche Ande- 
rung pm « dp/dt gebildef und 

f) die Compliance-Funktion C(p) mindestens aus p(t), 
p(nund R nach einem M dell berechnet wird. 

Die in-vivo bestimmte individuelle Compliance-Funktion 
C(p) kann zur Modellierung des Herzkreislaufsystems, 
insbes ndere zur kontinuieiiichen Berechnung des syste- 
mischen Blutflusses und des Herzzeitvolumens aus der 
Pulskurve p(t) verwendet warden. 
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Beschrcibung 

Die Erfindung betriffl ein Verfahren zur individuellen in-vivo Bestimmung der Compliance-Funklion C(p) = dV/dp 
dcs Gefaflsystems nach einer Herzkammer eincs Lebewesens aus dem Blutdruck p(t) und einem Referenz-Herzzeitvolu- 
5 menHZVref. 

Die Erfindung beirifft auch ein Verfahren zur laufenden Besummung des systemischen Blutflusses q(t), bei dem der 
Blutdruck p(t) in der Aorta oder in der Nahc der Aorta laufend bestimmt wird. 

Femer betrifft die Erfindung auch eine Vonichtung zur Durchfuhrung der Verfahren, bei der mindestens eine hamody- 
namische GroBe bestimmt und zur Anzeige, zum Ausdruck oder zur Weiterverarbeitung bereitgestellt wird, welche Vor- 

10 richtung mindestens einen Blutdmcksensor aufweist. 

Verfahren und Vorrichtungen der oben genannten Art sind bekannt. In der Vergangenheit wurden mehrere Verfahren 
mit dem Ziel entwickelt, das Herzzeitvolumen (HZV) aus dem arteriellen Blutdruck zu berechnen. Eincrseits werden 
Methoden vorgeschlagen, in denen das HZV aus wenigen charakterisuschen GrbBen wie etwa mittlerem arteriellen 
Druck (MAD), systolischem und diastolischem Druck (ADsys, ADdia), Systolendauer (ET) und Patientendaten (Alter, 

15 Geschlecht, Gewicht, GroBe) bestimmt wird [4, 6, 7], andererscits werden Algorithmen eingesetzt, in denen die vollstan- 
dige Kontur der pulsierenden Blutdruckkurve zur Berechnung des Herzzeitvolumens ausgenutzt wird [1,5, 20]. In den 
letztgenannten Methoden, die auch als Pulskonturanalyse bezeichnet werden, wurden dabei bisher zwei unterschiedliche 
Ansatze verfolgt. Zum einen wurde das HZV direkt aus dem arteriellen Blutdruck mit Hilfe einiger Korrekturfaktoren 
berechnet [19, 20], wohingegen in anderen Arbeiten [5, 25] aus dem Druck zusammen mit bestimmten Annahmen ein 

20 BlutfluB berechnet wird, von dem angenommen wird, daB er mit dem tatsachlichen BlutfluB in der Aorta ubereinstimmt 
und damit zur Bestimmung des Herzzeitvolumens verwendet werden kann. 

Das klassische Windkesselmodell, das erstmals von Hales [26] vorgeschlagen und von Frank [27] eingesetzt wurde, 
urn das Schlagvolumen (SV) und, zusammen mit der Herzfrequenz, das Herzzeitvolumen zu bestimmen, verwendet le- 
diglich den pcriphcrcn Widcrstand R und die Compliance C zur Modellierung des untersuchten Herzkreislauf systems. In 

25 diesem Modcll wird der artcricllc BluuluB durch q(t) beschrieben, der mit Hilfe des zu messenden Blutdruckes p(t) fur 
gegebencs C und R berechnet werden kann. Genauere Untersuchungen zeigen jedoch, dafl dieses cinfache Modell die 
physiologischcn Vcrhaltnissc nur un/.urcichend wiedergibt, so daB zahlreiche Modifikationen des urspriinglichen Mo- 
dells vorgeschlagen wurden; fiir cine Obcrsicht wird auf [24, 28] verwiesen. 

Die Genauigkcit dicscr Verfahren hangt dabei entscheidend davon ab, wie gut die Voraussetzungen, d. h. das verwen- 

30 dete Modcll, die Vcrhaltnissc im untersuchten Herzkreislauf system widcrspiegelt, so wird in [5] ein nicht lineares Wind- 
kesselmodell eingesetzt. dessen Parameter durch Alter und Geschlecht des Pauenten vorbesummt sind. Neuere Untersu- 
chungen [21 ] zeigen jedoch. daB das in |5| eingesetzte Modell nicht die korrekten physiologischen Randbedingungen 
wiedergibt; insbesonderc kann die Compliance (Dehnbarkeit) der GefaBe nicht stets durch den in [5] angegebenen Com- 
pliance-Druck-Zusammcnhang beschrieben werden. Fur diese Diskrepanz kommen mehrere Ursachen in Frage. Zum ei- 

35 nen wird in |5] nur cine Abhangigkcit dcs in-vitro bestimmten Aortenquerschnitts vom Blutdruck beriicksichtigt und 
eine Langenandcrung, wic in [22, 23] beschrieben, vernachlassigt; auch wird nur die Dichte des Blutes und nicht die 
stark hamatokritabhangigc Viskositat beriicksichtigt, ebenfalls bleibt die Compliance des peripheren Systems unberiick- 
sichtigt. Zum anderen kann das in [5] beschriebene Verfahren, abgesehen von Alter und Geschlecht, nicht die Compli- 
ance-Funkuon C(p) eincs cinzclnen Individuums verwenden. Gerade bei der Untersuchung pathophysiologischer Falle, 

40 z. B. ArterioskJerosc, kann jedoch nicht davon ausgegangen werden, daB sich C(p) dem Alter und Geschlecht entspre- 
chend verhalt, so daB dann das zugrundeliegende Modell die physiologischen Verhaltnisse nur unzureichend beschreibt 
[251. Letztlich wurde in [24] aufgezeigt, daB zu erwarten ist, daB ein modifiziertes Windkesselmodell die physiologi- 
schen Randbedingungen cxaktcr wiedergeben kann. 
Allen bisher diskuticrtcn Modellen ist jedoch gemein, daB die Modellparameter, nachdem sie einmal bestimmt wur- 

45 den nicht mehr vom Zustand dcs untersuchten Herzkreislaufsystems abhangen. Nahezu alle Parameter kdnnen sich je- 
doch mit der Zeit vcrandcrn, so kann sich beispielsweise der systemische Widerstand R auf grand emer Medikamentie- 
rung verandern. Anderc Parameter, so auch die Dehnbarkeit und Lange der Aorta, andern sich mit dem Druck so stark, 
daB sie bereits wahrend eincs Hcrzschlages als veranderlich betrachtet werden miissen. 

DaB Aortenimpedanz und Compliance nicht als konstant angenommen werden kimnen, wurde sowohl im Tierexpen- 

50 ment [22] als auch fiir Mcnschen gezeigt [29]. Ursache hierfiir sind primar die Dehnbarkeit, Langen- und Wumenande- 
rung der Aorta und aortennahen GefaBe. Der typische Verlauf des Aortenvolumens V in Abhangigkeit vom Druck wurde 
unter anderem in [30] beschrieben. Da die Compliance des Systems durch 

gegeben ist, muB aufgrund des begrenzten Volumens die Compliance fiir sehr hohe Driicke gegen Null streben und kann 
nicht konstant sein. Da die Volumenanderung durch Langen und Querschnittsanderung der GefaBe verursacht wird, an- 
dert sich ebenfalls die Aortenimpedanz, welche nach der Navier-Stokes-Gleichung zum einen von Querschnitt und 

60 Dichte p der Flussigkeit und zum anderen durch Lange, Viskositat t\ und Dichte des Blutes besummt wird. 

Druckabhangige Aortenimpedanz und Compliance wurden unter anderem in [5, 21] diskuuert und dort zur Untersu- 
chung nicht linearer Windkesselmodelle eingesetzt. In [5] wind dabei insbesondere davon ausgegangen, daB C(p) durch 
Alter und Geschlecht des Patienten festgelegt werden kdnnen. Auch bleibt die Impedanz Z(p) in diesem Ansatz unoe- 
riicksichtigt. Zudem folgt aus den in [21] erzielten Ergebnissen, daB das in [5] verwendete Modell teilweise im Wider- 

65 spruch zu den physiologischen Gegebenheiten steht. Ursache hierfiir ist zum einen die \forgabe derComphance und Aor- 
tenimpedanz. Ein solcher Ansatz ist ungeeignet, die fiir den untersuchten Patienten charakterisuschen Merkmale zu be- 
riicksichtigen. Auch kann das in [5] vorgestellte Verfahren nicht ohne Anderungen auf andere Spezies Ubertragen wer- 
den. Zum anderen wird in [5] nur der zuvor in-vitro untersuchte typische Aortendurchmesser und die Dichte des Blutes 



2 



DE 198 14 371 A 1 

beriicksichtigt. Der EinfluB von Langenanderungen der Aorta, sowie das dynamische \ferhalten der aortennahen und pe- 
ripheren GefaBe und die Viskositat des Bluies bieiben bei der Modellierung der in- vivo voriiegenden Bedingungen un- 
beriicksichtigt. 

Demzufolge existiert bisher kein Verfahren, das zur individuellen in- vivo Besummung des Compliance-Druck-Zu- 
sammenhangs die hier eingesetzten MeBgroBen verwendet. - ^ 

Diese Nachteile sollen durch das erfindungsgemaBe Verfahren ausgeschiossen werden, indem alle interessierenden 
Modellparameter aus Messungen am untersuchten physiologischen System, d. h. Mensch oder Tier, besummt werden. 
Hierzu ist insbesondere fortlaufend der aorten oder aortennahe Blutdruck p(t) und zu mindestens einem Zeitpunkt ein 
Referenz-Herzzeitvolumen (HZVref) zu messen. Mit Hilfe dieser GroBen kdnnen alle Parameter festgelegt werden und 
anschlieBend zur Untersuchung der Hamodynamik eingesetzt werden. i 10 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur individuellen in-vivo Bestimmung der Compliance- 
Funktion eines Lebewesens zu schaffen, die die physiologischen Verhaltnisse moglichst getreu wiedergibt. 

Der Erfindung liegt des weiteren die Aufgabe zugrunde, ein Verfahren zur laufenden Bestimmung des systemischen 
Blutflusses eines Pauenten zu schaffen, das gering invasiv ist und den tatsachlichen Blutflufl zu jedem Zeitpunkt mog- 
lichst exakt beschreibt. is 

Ferner liegt der Erfindung die Aufgabe zugrunde, eine Varrichtung zur Durchfuhrung solcher Verfahren zu schaffen, 
mil deren Hilfe eine zuverlassige individuelle Bestimmung der Comphance-Funktion und eine laufende zuverlassige 
Uberwachung des systemischen Blutflusses eines Pauenten unter minimal invasiven Bedingungen moglich ist. 

Die Losung dieser Aufgaben ist in den unabhangigen Patentanspriichen angegeben. \farteilhafte Ausfuhrungsformen 
sind den abhangigen Patentanspruchen zu entnehmen. 20 

GemaB der Erfindung wird ein sogenanntes Windkesselmodell verwendet, dessen Parameter unter Zuhilfenahme einer 
Referenz-Messung in-vivo identifiziert werden konnen. AnschlieBend werden damit der systemische FluB und weitere 
hamodynamische Parameter besummt. Ein derart angepaBtes und modifiziertes Windkesselmodell beschreibt das Herz- 
kreislaufsystem des untersuchten Individiums genauer und kann somit zur ebenfalls genaueren Berechnung des systemi- 
schen Flusses und daraus abgeleiteter hamodynamischer Parameter verwendet werden. Auch kann das Verfahren unmit- 25 
telbar auf andere Spezies ubertragen werden, ohne hierfur zuvor hypothetische Charakteristika zu bestimmen. Der Mehr- 
aufwand fur die neu entwickelte Methode besteht darin, daB zur Berechnung einer individuellen Compliance-Funkuon 
neben der kontinuierlichen Blutdruckmessung mit einem anderen Verfahren, z. B. arterieller Thermodilution, das Herz- 
zeitvolumen fur mindestens einen Zeitpunkt bestimmt werden muB. 

Die Erfindung wird im folgenden an Hand eines in den Figuren schemausch dargestellten Ausfuhrungsbeispiels naher 30 
erlautert. Es zeigt: 

Fig. 1 Eine bevorzugte elektrische Modellschaltung zur Simulation des untersuchten Herzkreislaufsystems. 
Fig. 2 Ein FluBdiagramm zur Berechnung der Aortenimpedanz Z(p), der Compliance-Funkuon C(p) und des Blutflus- 
ses q(t). 

Fig. 3a Ein Schaubild zur Erlauterung der zeithchen Abhangigkeit des Blutdrucks p(t), wobei to den Zeitpunkt des 35 
Offhens der Aortenklappe, t^ den Zeitpunkt des SchlieBens der Aortenklappe und t D den Zeitpunkt des Endes der Dia- 
stole kennzeichnen. 

Fig. 3b Ein Schaubild zur Erlauterung der zeitlichen Abhangigkeit des Blutdmckes p(t) und des resultierenden Blut- 
flusses q(t) 

Fig. 4 Eine typische Compliance-Funkuon C(p) der Aorta eines Menschen. 40 
Fig. 5 Ein Blockschaltbild einer Vorrichtung gemafl der Erfindung. 

Fig. 1 zeigt ein nicht lineares und modifiziertes Windkesselmodell, das erfindungsgemaB bevorzugt eingesetzt wird, 
worm die Aortenimpedanz-Funkuonen Z(p) und Zo(p), die Compliance-Funktion C(p) und der systemische Widerstand 
R berucksichugt werden. 

Der Widerstand R in Fig. 1 reprasentiert den zeitlich variablen peripheren Stromungswiderstand des Korpers. Zo(p) 45 
und Z(p) sind nicht linearc, von dem Druck p(t) abhangige Impedanzen, die zusammen mit der nicht linearen, druckab- 
hangigen Kapazitat C(p) das Verhalten der Aorta und der aortanahen BlutgefBBe simulieren sollen. 

Fiir das in Fig. 1 skizzierte Modell resultiert mit der Fouriertransfonnierten 



p(o)) » J p(t)exp(-icot)dt 

und der analog zu berechnenden Funkuon q (cu) fur die Aortenimpedanz 
o^co p(o>) n L 



50 
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so daB fur R » Z folgt 



Z= lim 



P((0) 
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(2) 
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Fur die Compliance C gilt dariiber hinaus 

cm qft)-P(t)/R 

° lP; ~ p(t)-Z(p).-(q(t)-p(t)/R) , (3) 

5 .. 

darin sind p(t) = dp(t)/dt und q(t) = dq(t)/dt die zeitlichen Anderungen des Druckes bzw. Flusses. Gleichungen (1) bis (3) 
zeigen, dafi C und Z fiir das Modell aus Fig. 1 berechnet werden k6nnen, wenn der systemische FluB q(t), der Blutdruck 
p(t) und der systemische Widerstand R bekannt sind. Das Modell wird jedoch das untersuchte Herzkreislaufsystem nur 
io so gut beschreiben, wie dieses im Rahmen des vorgestellten Ansatzes moglich ist. 
Einen Uberblick uber das entwickelte Verfahren zeigt Fig. 2. 

i) Zunachst wird aus dem Druck p(t) der mittlere arterielle Blutdruck MAD und die Herzrate HR bestimmt. 

ii) Zusammen mit dem Referenz-Herzzeitvolumen HZVref, das vorzugsweise durch arterielle Thermodilution be- 
15 stimmt wurde, und fiir das gilt 

HZVref = HR • Jq(t)dt (5) 

wird der systemische Widerstand nach R = (MAD-ZVD)/HZVref berechnet. 
20 Darin ist ZVD der zentralvenose Druck, der, falls er unbekannt ist, durch einen geeigneten konstanten Druck, z. B. 

0 mm Hg, approximiert werden kann. 

iii) Als nachstes ist ein geeignet zu wahlender FluB q(t) festzulegen, der als Startfunktion der folgenden Iteration 
diem und den physiologischen Randbedingungen genugt. Der BlutfluB q(t) beschreibt denjenigen FluB, der unmit- 
telbar vom Linksherz in die Aorta hineinstromt. Deshalb wird fur q(t) gefprdert, daB die Nebenbedingungen 

25 

HZVref 

q(t 0 )=q(ts) - 0 und Jq(t)dt = — j^p = SV ( 6 ) 

fur alle Herzschlage erfullt ist, die wahrend der Messung des Referenz-Herzzeitvolumens HZVref registriert wur- 
den. Der Zeitpunkt !<, ist darin der Systolenbeginn und t, das Systolenende. Das Ende der Diastole wird im folgen- 
den mit t D bezeichnet. Geeignete Fliisse q(t) waren q(t) = 0 fin- 1* < t < t D und fur to < t < t, etwa 

q(t ) g *' H ™*\ sin (pSp (t-to ) ) . (7) 
HK ' 2-HR-(ts-to) l s' to 

oder bevorzugt 

.o q(t) = K[p(t)-(p(t 0 ) + (t-t 0 ) P(ts t ^ tQ) ) l . (8) 

bei letzterem ist dabei K aus der Bedingung Jq(t)dt = HZVref/HR zu bestimmen. Neben den in Gleichung (7) und 
(8) angegebenen Flussen sind jedoch auch andere Anfangsbedingungen, wie etwa ein konstanter oder parabelartiger 
45 Flup, denkbar. 

iv) Setze Z(p) = 0 und filhre die Hilfsvariablen qai t (t) und ein, die mit q^t) = q(t) und E^-co miualisiert wer- 
den. 

v) Berechne aus dem aorten oder aortennah gemessenen Blutdruck p(t), dem FluB q(t) und deren zeitlicher Ablei- 
tungen die Compliance-Funkuon nach 
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50 



c, D) _ q(t)-P(t)/R 

MP ' p(t)-Z(p)-(q(t)-p(t)/R) . (9) 

55 vi) AnschlieBend wird die Inverse der Compliance-Funktion durch ein Polynom geeigneter Ordnung in p approxi- 
miert, d. h. in 

— = Ta k P k ( 10) 
C(p) ^ akP 
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werden die Entwicklungskoeffizienten ax so bestimmt, daB die Identity 

q(t) . _P0) + J_ -lp(t) - Z(p) • (q(t) - p(t)/R)] 

QU R t2ic«kP 1 (11) 

optimal erfullt wird. Als geeignetes Kriterium hierzu kann der quadrausche Fehler 
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E = <(q(t)-p(t)/R-[p(t)-Z(p) • (^tHW/RMZkOkp 11 ]) 2 ) 

minimiert werden, wobei p(t) und q(t) zu alien Zeitcn oder auch nur aus bevorzugten Zeitintervallen verwendet wer- 
den konnen. Die Notation ( • ) bedeutet dabei hier und im folgenden die Berechnung des Mittel wertes. 

vii) Falls E < E^,, dann setze q^l) = q(t) und E^i = E und fahre mil Schritt viii) fort, anderenfalls gehe zu Punkt x). 

viii) Berechne Z(p). Zum einen kann hierzu so vorgegangen werden, daB Z(p) nach Gleichung (1) Oder (2) zu be- 
stimmen ist. Dabei wird angenommen, daB diese Gleichungen, die zunachst nur fur das Modell aus Abb. lb gel ten, 
in dem die Parameter nicht druckabhajigig sind, fur geniigend kurze Zeitintervalle At auch fur den hier vorgestellten 
nicht lineare Ansatz entsprechend Fig. 1 zutreffen. FUr diesen kann dann mit 

+At 

Z(p- Jp(t)dt/2At) 

-Al 

die Impedanz-Funktion nach Gleichung (2) oder gleichwertig durch 15 
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/ t+At l+At t+At t+At 

Z(P) =l/( / P 2 (t)dt- ( J p(t)dt) 2 )/( f q 2 (t)dt - ( / q(t)dt) 2 ) 

I t-At t-At t-At t-At (12) 20 

auch direkt aus dem zeilabhangigen Blutdmck p(t) und dem BlutfluB q(t) ohne vorherige Fouriertransfonnation er- 
mittelt werden. 

Andererseits kann Z(p) jedoch noch einfacher durch 25 

(13) 

30 

berechnet werden, wic folgcndc Uberlcgung zeigt. Aortendurchmesser d und -Lange I nehmen mit zunehmendem 
Druck zu, so daB in niedrigstcr Nahcrung d a I angesetzt wird. Nach dem Hagen-Poiseuilleschen Gesetz resultiert 
damit Z(p) a T|/V, worin r\ die Viskositat des Blutes und V das Aortenvolumen bezeichnet. Mit C(p) = dV/dp ergibt 
sich damit C(p) a d(l/Z)/dp, woraus Gleichung (13) unmitteibar folgt. Die darin enthaltene Proportionalitatskon- 
stante A kann z. B. dadurch bcsiimmt werden, daB fur mindestens einen Druck p die Funktion Z(p) nach Gleichung 35 
(I2)bestimmt wird. 

ix) Falls der Fehler E geniigend klein ist, dann wird die Identification der Modellparameter hier abgebrochen. An- 
derenfalls wird mit x) fortgefahren. 

x) AbschlieBend ist der angenommene BlutfluB derart zu Variieren, daB das Schlagvolumen weiterhin dem Schlag- 
volumen SV = HZVref/HR cntspricht, das aus dem Referenz-Herzzeitvolumen folgt. Da zu diesem Zeitpunkt q^l) 40 
stets den bis dahin optimalen FluB beschreibt, wird q(t) = q a j l (t) + 5q(t) mit J5q(t)dt = 0 gesetzt. 

xi) Fahre mit Schrift v) fort. 

Der in i)-xi) angegebene Algorithmus beschreibt das bevorzugte Verfahren, in dem aus dem Referenz-Herzzeitvolu- 
men HZVref und der laufend gemessenen arteriellen Pulskurve p(t) alle Ubrigen GriJBen bestimmt werden. Hierdurch 45 
wird gewahrleistet, daB die Compliance-Funktion und Aortenimrjedanz-Funktion so bestimmt werden, wie es die effek- 
uv in dem untersuchten Herzkreislaufsystem wirkenden Zusammenhange erfordern. Insbesondere benicksichtigt C(p) 
dadurch nicht nur die Anderung des Aortenquerschnitts, sondern die tatsachliche \folumenanderung von Aorta und peri- 
pheren GefaBen; ebenso wird durch Z(p) eine Langenanderung der Aorta und die Dichte und Viskositat des Blutes be- 
schrieben. Mit Hilfe der in Schritt vi) eingefuhrten Beschreibung, die jeden physiologisch moglichen Compliance- 50 
Druck-Zusammenhang wiedergeben kann, wird zudem die Extrapolation von C(p) und Z(p) ermdglicht, so dafi die Funk- 
tionen auch uber die wahrend der Kalibrierung, d. h. iiber die wahrend der Messung des Referenz-Herzzeitvolumens be- 
obachteten Driicke hinaus angewendet werden konnen. 

Der angegebene Algorithmus beschreibt das bevorzugte Verfahren. Andere Verfahren konnen leicht aus diesem ge- 
wonnen werden und werden durch die aufgefuhrten Anspruche mit abgedeckt. Insbesondere kann in vi) statt des Kehr- 55 
wertes der Compliance-Funktion auch C(p) durch eine endliche Taylorreihe, d. h. durch ein Polynom, beschrieben wer- 
den. Auch ist denkbar, dafi die Schritte viii) und ix) miteinander vertauscht werden oder, daB in Schritt vi) das Kriterium 
fttr die Opurnierung gcandert wird. So konnte statt der Minimierung des quadratischen Fehlers auch eine Maximierung 
des Erwartungswertes erfolgen. 

Urn das Verfahren zu beschleunigen, kann insbesondere die Wahl des anfanglichen Blutflusses q(t) in Schritt iii) derart 60 
gestaltet werden, dafi C(p) zunachst nur aus den Diastolen, d. h. fur t, < t £ to bestimmt wird. Nach Gleichung (9) gilt 
dort fur Z = 0 dann C(p) = -p(t)/R • p(t), so daB hiermit fur alle p(t) < p(ts) die Compliance-Funktion dazu ausgenutzt 
werden kann, den FluB wahrend der Systole nach einem orthogonalen Funktionensystem, z. B. 

q(0=I k q k sin(kMi.(t-t o )/(t s -t o )) (14) 65 
zu entwickeln. Die Koeffizienten qk werden dabei durch Minimierung des Ausdrucks 
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(((q(l)-p(t)/R-p(l)/C(p)) 2 ) (15) 

bestimmt, wobei der obige Fehler aus alien p(t) < p(ts) zu berechnen ist, einschlieBlich derjenigen, die zum Beginn der 
Systole auftreten. Mit dieser Wahl des Flusses geniigen bei einem ausreichend guiartigen Herzkreislaufsystem bereits 
5 eine oder wenige Iterauonen, um nach dem beschriebenen Algorithmus die Modellparameter zu berechnen. 

Eine weitere Beschleunigung des Verfahrens kann fur den Fall erfolgen, daB in Schritt vi) nur Druckwerte aus den Be- 
reichen der Diastole verwendet werden und 1/C(p) durch ein Polynom zweiten Grades approximiert wird. Wahrend der 
ersten Iteration resultiert fur diese Anwendung 

io MAD 2 * HZVref 1 

(p(t)> 3 • MAD • p - 3 'MAD 2 - p 2 (16) 



Nachdem die Parameter des Modells aus Fig. 1 aus dem aorten oder aortennahen Blutdruck p(t) und dem Referenz- 
15 Herzzeitvolumen HZVref bestimmt wurden t kann das Herzzeitvolumen anschlieBend kontinuierlich alleine aus dem 
Blutdruck von Schlag zu Schlag berechnet werden. Hierzu ist der Beginn der Herzschlage, die im folgenden wiederum 
mit ^ bezeichnet werden, zu berechnen und Gleichung (11) fiir die Anfangsbedingung q(to) = qfa) = 0 zu berechnen. 
Eine solche Berechnung kann numerisch erfolgen, z. B. nach Runge-Kutta, Euler oder einem anderen allgemein bekann- 
ten Algorithmus. Liegt die erforderliche Blutdruckmessung dazu nur in diskreten Zeitschritten vor, so kann zwischen 
20 diesen bei Bedarf linearisicrt werden. Fiir das so nur mit Hilfe des Blutdrucks bestimmte Schlagvolumen gilt mit y = 1 

SV = yJq(t)dt (17) 

so daB fiir das kontinuierlich bestimmte Herzzeitvolumen HZV folgt 

25 

HZV = HR SV = HR • y ■ Jq(i)dt (18) 

Darin bezeichnet HR die Hcrzratc, die ebenfalls aus der Pulskurve p(l) zu berechnen ist und die Zahl der Schlage pro 
Minute angibt. Die in Gleichung (16) bzw. (17) auftretende Integrauon kann dabei uber den gesamten Herzschlag erfol- 

30 gen oder auch nur Ubcr den Zcitraum der Systole, da wahrend der Diastole gilt q(t) 0. Wird das Schlagvolumen SV, 
und damit auch HZV aus dem BlutlluB q(t) wahrend des gesamten Herzschlages berechnet, so muB nicht das Systolen- 
eridc bestimmt werden. Hierzu ware (siehe z. B. [31 ]) cntweder eine genaue Analyse der Druckkurve notwendig, um aus 
p(t) die Lage des sogenanntcn Dichrotic Notch und damit das Systolenende zu bestimmen, oder weitere MeBvorrichtun- 
gen, wie etwa ein EKG, erfordcrlich. Die Integrauon uber den gesamten Zeitraum ist demzufolge robuster und weniger 

35 aufwendig als dicjenigen Vcrfahren, die nur einen bestimmten Zeitraum des Herzschlages auswexten. Wird zudem das 
kontinuierlich bestimmte Herzzeitvolumen HZV auch aus denjenigen Blutdruckmessungen berechnet, die zusammcn 
mit dem Referenz-Herzzeitvolumen HZVref zur Modellidentifikation nach dem oben beschriebenen Verfahren verwen- 
det wurden. so kann die Genauigkeit des Verfahrens zur kontinuierlichen HZV-Berechnung noch dadurch gesteigert wer- 
den, daB HZV = HZVref gefordert wird und damit der Kalibrierfaktor y nach 

HZVref (19) 
Y HR • fq(t)dt 



45 bestimmt wird. . . . . . 

Um das beschriebene Verfahren einzusetzen, ist eine \forrichtung erforderlich, deren pnnapieller Aufbau in *ig, 4 
dargestellt ist. Die darin gestrichelt dargestellten Komponenten sind optional und konnen in einer Minimalkonfigurauon 
der Vorrichtung zumindest teilweise entfallen. Eine solche Vorrichtung besteht mindestens aus einer Auswerteeinheit, 
gewbhnlich einer Recheneinheit, in der das Verfahren zur Bestimmung der individuellen CompUance-Funktion C(p) al- 
so leine oder zusammen mil anderen Methoden eingesetzt wird; insbesondere kann das Verfahren zur kontinuierlichen Be- 
rechnung des Herzzeitvolumens in der selben Vorrichtung verwendet werden. Zudem sind erforderlich ein Sensor zur 
Messung des aorten oder aortennahen Blutdrucks p(t) und eine Signalaufbereitung und Signalwandlung, ein Programm- 
und ein Datenspeicher, sowie eine Vorrichtung zur Vorgabe des Referenz-Herzzeitvolumens HZVref. Wird das HZVref 
uber eine arterielle Thermodilution bestimmt, so besteht diese Einheit mindestens aus einem Bluttemperatursensor und 
55 einem Sensor zur Messung der Temperatur des Injektats, welches von diesem Verfahren verwendet wird, siehe [8]. Da 
HZVref jedoch auch auf anderem Wege vorgegeben werden kann, kann dieses Modul auch anders ausfallen oder eine 
Eingabe liber Tastatur erfolgen, die zudem zur Benutzerfuhrung in der Vorrichtung verwendet werden kann. Zudem be- 
steht mindestens eine der MGgkchkeiten, die Ergebnisse der Auswertung anzuzeigen, auszudrucken oder auf emem Mas- 
senspeicher (nicht dargestellt) zu speichem. 
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Patentanspriiche 

1. Verfahren zur individuellen in-vivo Bestimmung der Compliance-Funktion C(p) = dV/dp des GefaBsystems 
nach einer Herzkammer eines Lebewesens aus dem Blutdruck p(t) und einem Referenz-Herzzeitvolumen HZVref, 20 
dadurch gekennzeichnet, daB 

a) der Druck p(t) in der Aorta oder in der Nahe der Aorta laufend bestimmt wird; 

b) ein Referenzwert fur das Herzzeitvolumen HZVref verwendet wird; 

c) aus dem Blutdruck p(t) der mittlere Blutdruck MAD berechnet wird; 

d) der systemische Widerstand R des Korpers berechnet wird als 25 

MAD - ZVD 
H " HZVref 

wobei ZVD ein beliebiger zentralvenoser Druck ist, der ermittelt oder geschatzt wird; 30 

e) aus dem Blutdruck mindestens die erste zeithche Anderung p(t) = dp/dt gebildet wird; und 

0 die Compliance-Funktion C(p) mindestens aus p(t), p(t) und R nach einem nicht linearen Modell berechnet 
wird. 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dafl das Referenz-Herzzeitvolumen HZVref mittels einer 
Dilutionsmethodc bestimmt wird. 35 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB als Diluuonsmethode eine Thermodiluuonsmethode 
eingesetzt wird. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daB eine arterielle Thermodiluuonsmethode eingesetzt 
wird. 

5. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB als Diluuonsmethode eine Farbstoffdiluuonsmethode 40 
eingesetzt wird. 

6. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB fur die Berechnung der Compliance-Funktion C(p) 
nur Werte von p(t) verwendet werden, fur die gilt: 

p(0 < p(0, 45 
wobei t* der Zeitpunkt des Schliefiens der Aortenklappe ist. 

7. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB fur die Berechnung der Compliance-Funktion C(p) 
nur Blutdruckwerte aus der Diastole verwendet werden. 

8. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB fur die Berechnung der Compliance-Funktion C(p) 50 
nur Blutdruckwerte aus der Systole verwendet werden. 

9. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB aufgrund des Druckes p(t) 
und der ersten zeitlichen Ableitung dp/dt ein BlutfluB q(t) bestimmt wird und die Compliance-Funktion nach 

c( P )= qw-rttyn 

p(t)-Z(p)-(q(t)-p(t)/R) 
fur beliebige Impedanzenfunktionen Z(p) und beliebige Zeiten t derart berechnet wird, daB 

q(t) = + C(p)[p(t) - 2(p) • (q(t) - p(t)/R)] 60 

bestmdglich erfullt wird. 

10. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB der Kehrwert der Compli- 
ance-Funktion, d. h. 1/C(p), durch ein Polynom endlicher Ordnung approximiert wird und dieses Polynom zur Ex- 65 
trapolation von C(p) uber den wahrend der Bestimmung des Referenz-Herzzeitv lumens erfaBten Druckbereich 
hinaus verwendet wird. 

11. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB das Minimum der Funktion 
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£ j s ( q(t) . m efiLna jqto-pw/R) ) 2 nt 



wo •• **°*p 

bestimmt wird und dafl die individuelle Compliance-Funktion C(p) danach berechnet wird als 

C(P) = - 1 k 
2k«kP 

12. Verfahren nach Anspruch 7 und 10, dadurch gekennzeichnet, dafl der Kehrwert der Compliance-Funktion C(p) 
durch ein Polynom zweiter Ordnung beschrieben wird und C(p) durch folgende Funktion approximiert wird 

. MAD 2 - HZVref 1 

<P(t)> 3 • MAD -p - 3 -MAD 2 - p 2 

1 3. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB die Compliance-Funktion 
C(p) durch ein Polynom endlicher Ordnung approximiert wird und dieses Polynom zur Extrapolation von C(p) iiber 
den wahrend der Bestimmung des Referenz-Herzzeitvolumens erfaflten Druckbereich hinaus verwendet wird. 

14. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl das Minimum der Funktion 



2/"(q(t) " ' Ek PkP k ] lP(t) " Z(p) • (q(t) - p(t)/R)]) 2 dt 



t-to 

bestimmt wird und dafl die individuelle Compliance-Funktion C(p) danach berechnet wird als 
C(p) = I k ft t p k . 

15. Verfahren nach Anspruch 7 und 9, dadurch gekennzeichnet, dafl die fur p(t) < p(t s ) berechnete Compliance- 
Funktion C(p) zur EntwickJung des Blutflusses q(t) nach einem vollstandigen Funktionensystem verwendet wird; 
insbesondere wird q(t) in Form einer Fourierreihe durch folgende Gleichung beschrieben 

q(t)-£q k sin (k-ir-gj;) , 
k 

wobei die Koeffizienten q^ durch Minimierung des quadratischen Fehlers bestimmt werden und die Groflen t<, und ts 
die Zeitpunkte des Offnens und Schlieflens der Aortenklappe bezeichnen. 

16. Verfahren nach einem der vorhergehenden Anspniche, dadurch gekennzeichnet, dafl der angenommene Blut- 
flufl q(t) derart variiert wird, dafl der quadrausche Fehler minimal ist. 

17. Verfahren nach einem der vorhergehenden Anspniche, dadurch gekennzeichnet, dafl der Aortenimpedanz- 
Druck-Zusammenhang durch 



bestimmt wird, wobei A eine Proportionalitatskonstante ist. 

18. Verfahren nach einem der vorhergehenden Anspniche, dadurch gekennzeichnet, dafl eine nicht lineare Aorten- 
impedanz-Funktion mil Hilfe der Fouriertransformierten von Blutdruck p(<o) und postuliertem Blutfiufl q((u) nach 
1/Z =q(a>- > oo)/p(a>- > »)-l/R oder Z(p) = p(u>- > »)/q(a>- > oo) bestimmt wird. 

19. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl aufgrund des Dmckes p(t) und der ersten zeitlichen 
Ableitung dp/dt ein Blutfiufl q(t) bestimmt wird und die Impedanz-Funktion 

+At 

Z(p- Jp(t)dt/2At) 

-Al 

nach 

/ t+Ai uai t*At t+At 

Z(P) A ( J P 2 (t)dt- ( / p(t)dt) 2 )/( / q 2 (t)dt - ( J q(t)dt) 2 ) 

V [-At t-At t-At t-At 
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berechnet wird. 

20. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, dafi die Aortenimpedanz- 
Funktion Z(p) durch ein Polynom endlicher Ordnung approximiert ward und dieses Polynom zur Extrapolation von 
Z(p) liber den wahrend der Kalibricrung erfafiten Druckbereich hinaus verwendet wird. 

21 . Verfahren zur laufenden Bestimmung des systemischen Blutflusses q(t), bei dem der Blutdruck p(t) in der Aorta 5 
oder in der Nahe der Aorta laufend bestimmt wird, dadurch gekennzeichnet, daB mindestens eine in- vivo besummte 
individuelle Compliance-Funkuon C(p) zur Modellierung des Herzkreislaufsystems verwendet wird. 

22. Verfahren nach Anspruch 21, dadurch gekennzeichnet, daB der systemische Blutflufl q(t) besUmmt wird mit: 

q(t)- -^ + C(p)[p(t)-Z(p) -(q(l)-p(t)/R)] 10 

R . 

23. Verfahren nach Anspruch 21, dadurch gekennzeichnet, daB die nach einem der Anspriiche 1 bis 20 ermittelte 
Compliance-Funkuon C(p) verwendet wird. 

24. Verfahren nach Anspruch 21 , dadurch gekennzeichnet, daB das Schlagvolumen SV durch Integration des Flus- 15 
ses uber einen geeigneten Zeitraum entsprechend 



SV = Jq(t)dt 

berechnet wird, wobei der geeignete Zeitraum insbesondere der Dauer des Herzschlages oder der Systolendauer 20 
wahrend des Herzschlages entsprechen kann. 

25. Verfahren nach Anspruch 1 und 21, dadurch gekennzeichnet, daB durch Vergleich des kontinuierlichen Blut- 
fiusses q(t) mit einem Referenz-Herzzeitvolumen HZVref das Schlagvolumen SV durch 

r HZVref 25 

SV = Yjc|(tdt mit y = : 

HR-Jq(t)dt 



berechnet wird. 

26. Verfahren nach einem der Anspriiche 24 bis 25, dadurch gekennzeichnet, dafi die Schlagvolumenvariauon nach 30 



sw =Y<sv 2 >-(sv) 2 

berechnet wird und alleine oder mit weiteren Parametem, z. B. mittlerem Blutdruck MAD, systolischem Druck AD- 35 
sys» diastolischem Druck AD DJA und Herzrate HR, zur Korrektur des Schlagvolumens verwendet wird. 

27. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB zusammen mit der Herz- 
rate HR aus dem Schlagvolumen SV konunuierlich das Herzzeitvolumen nach HZV = HR • SV berechnet wird. 

28. Verfahren nach Anspruch 27, dadurch gekennzeichnet, daB aus der Blutdruckkurve konunuierlich der mittlere 
Druck MAD bestimmt wird und damit der systemische Widerstand konunuierlich nach 40 

0 MAD - ZVD 



fur einen beliebigen zentralvenosen Druck ZVD berechnet wird, der gemessen oder geschatzt wurde. 45 
29. Verfahren nach Anspruch 27, dadurch gekennzeichnet, daB aus der Blutdruckkurve kontinuierlich der mittlere 
Druck MAD bestimmt wird und damit die Compliance-Funktion kontinuierlich nach 

C (p)= MAD 2 * HZV 1 

<P(t)> 3 • MAD -p - 3 -MAD 2 - p 2 50 



berechnet wird. 

30. Verfahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, dafi der Druck p(t) aus einem 
invasiven aortennahen Blutdruck, z. B. A.femoralis, A.brachialis, A.axillaris oder A jadialis bestimmt wird. 55 

31. Verfahren nach einem der vorhergehenden Anspriiche, in dem eine nicht invasive Druckmessung zusammen 
mit geeigneten Filtern, insbesondere Notch-Filter, zur Approximation des Druckes p(t) verwendet wird. 

32. Vorrichtung zur Durchfuhrung der Verfahren nach einem der vorhergehenden Anspriiche, bei der mindestens 
eine hamodynamische Grdfie bestimmt und zur Anzeige, zum Ausdruck oder zur Weiterverarbeitung bereitgestellt 
wird, welche Vorrichtung mindestens einen Blutdrucksensor aufweist, dadurch gekennzeichnet, dafi die Munch- 60 
tung eine Einrichtung zur in-vivo Bestimmung der Compliance-Funktion aufweist 

33. Vorrichtung nach Anspruch 32, dadurch gekennzeichnet, dafi die Einrichtung zum Bestimmen der Compliance- 
Funktion Mittel zur Bestimmung eines Referenz-Herzzeitvolumens oder Mittel zur manuellen Eingabe desselben 
umfafit. 

34. Vorrichtung nach Anspruch 32 oder 33, dadurch gekennzeichnet, dafi sie eine Auswerteeinheit aufweist, die 65 
mindestens den systemischen FluB q(t) aufgrund p(t) und der ersten zeitlichen Ableitung dp/dt berechnet. 

35. Vorrichtung nach Anspruch 34, dadurch gekennzeichnet, dafi der systemische Flufl q(t) berechnet wird mit 
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-^- ) + C(p)[p(t)-Z(p) -(q(t)-p(t)/R)] 
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